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We present results from a search for the flavor-changing neutral current decays B Kt^ and 
B -> K*£+£', where £+£~ is either an e^e or fj,'^ n pair. The data sample comprises 22.7 x 10 
T{4S) BB decays collected with the BABm. detector at the PEP-II B Factory. We obtain the 
90% C.L. upper limits B{B K£+£-) < 0.50 x 10"^ and B{B K'l+£-) < 2.9 x 10"^ close 
to Standard Model predictions for these branching fractions. We have also obtained limits on the 
lepton-family-violating decays B Ke^jj^ and B K*e^fi^. 

PACS numbers: 13.25.Hw, 13.20.He 



The flavor-changing neutral current decays B with branching fractions predicted to be of order 10^^ — 

K£+t~ and B K* {^'•)2)£+ 1~ , where i'^ is a charged 10~^ (l| ||. The dominant contributions arise at the 
lepton, are highly suppressed in the Standard Model, one-loop level and are known as electroweak penguins. 



4 



Besides probing Standard Model loop effects, these rare 
decays are important because their rates and kinematic 
distributions are sensitive to new, heavy particles — such 
as those predicted by supersymmetric models — that can 
appear virtually in the loop ||. 

The Standard Model predictions for B -> K^*H+l- 
include three main contributions: the electromagnetic 
(EM) penguin, the Z penguin, and the W^W~ box dia- 
gram. Evidence for the EM penguin amplitude has been 
obtained from the observation of -B — > K*^ and inclu- 
sive B — > Xs7, where Xg is any hadronic system with 
strangeness [|[ ||. 

Calculations of decay rates for B based 
on the Standard Model have significant uncertainties due 
to strong interactions. For example, Ali et al. predict 
B{B Ke+£-) = i0.57toll) x lO^^ for both e+e" and 
At+A*- final states, S(S ^ K*e+e-) ^ (2.3tj]J) x 10-^ 
and B{B K*fj.+fj.-) = (1.9tg:^) x 10"®- The contri- 
bution of the EM penguin amplitude to B ^ K*£^£^ is 
particularly strong at low values of m^+£- , giving a larger 
rate for B K*e+e^ than for B K*fi~^^^. 

We search for the following decays: 5+ K^i^i^, 
B° K'^£+£-, B+ K*+e+£-, and B" K*^l+t-, 
where K*^ K+tt- , K*+ K°n+, tt+tt- , and £ 

is either an e or ^. We also search for the lepton-family- 
violating decays B — > K'^*^e^^r^. Throughout this pa- 
per, charge-conjugate modes are implied. 

The data used in the analysis were collected with the 
BABAR detector at the PEP-II storage ring at the Stan- 
ford Linear Accelerator Center during 1999-2000. We 
analyzed a 20.7 fb~^ data sample taken on the T{AS) 
resonance consisting of (22.7 ± 0.4) x 10^ T(4S') BB 
events. 

This search relies primarily on the charged-particle 
tracking and particle-identification capabilities of the 
BABAR detector [|j . Charged particle tracking is provided 
by a five-layer silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH). The DIRC, a Cherenkov ring- 
imaging particle-identification system, is used for charged 
hadron identification. Electrons are identified using 
the electromagnetic calorimeter (EMC), which comprises 
6580 thallium-doped Csl crystals. These systems are 
mounted inside a 1.5 T solenoidal superconducting mag- 
net. Muons are identified in the instrumented flux return 
(IFR), in which resistive plate chambers are interleaved 
with the iron plates of the magnet flux return. 

We extract the signal using the kinematic variables 

'^ES = v/£^b'-(E,P:)' and = ^mi + pf - 
E^, where E^ is the beam energy in the e+e rest (cm.) 
frame, p* is the cm. momentum of daughter particle i in 
the B meson candidate, and is the mass of particle i. 
For signal events, tties peaks at the B meson mass with a 
resolution of about 2.5 MeV/c^ and Ai? peaks near zero, 
indicating that the candidate system of particles has total 
energy consistent with the beam energy in the cm. frame. 



To prevent bias in the analysis, we optimized the event- 
selection criteria using Monte Carlo samples: we did not 
look at the data in the signal region or in the sidebands 
that were used to measure the background until these 
criteria were fixed. Signal efficiencies were determined 
using the Ali et al. model |l]. 

We select events that have at least four charged tracks, 
the ratio i?2 of the second and zeroth Fox- Wolfram mo- 
ments 1^ less than 0.5, and two oppositely charged lep- 
tons with momentum p > 0.5 (1.0) GeV/c for e (//) can- 
didates. Electron-positron pairs consistent with photon 
conversions in the detector material are vetoed. We re- 
quire charged kaon candidates to be identified as kaons 
and the charged pion in K* — > Ktt not to be identified as 
a kaon. For B — > K*i^i^ , we require the mass of the K* 
candidate to be within 75 MeV/c^ of the mean isT* (892) 
mass. Kg candidates are reconstructed from two oppo- 
sitely charged tracks that form a good vertex displaced 
from the primary vertex by at least 1 mm. 

The decays B ^ J /il>{^ 1+ £-)K^*^ and 
B i!{2S){^ l+er)K(*'^ have identical topologies 
to signal events. These backgrounds are suppressed by 
applying a veto in the ls.E vs. rui+i- plane (Fig. |^). 
This veto removes charmonium events not only with re- 
constructed rrii+g- values near the nominal charmonium 
masses, but also events that lie further away in m£+£- 
due to photon radiation (more pronounced in electron 
channels) or track mismeasurement. Removing all of 
these events simplifies the description of the background 
shape. Charmonium events can, however, pass this 
veto if one of the leptons (typically a muon) and the 
kaon are misidentified as each other. If reassignment 
of particle types results in a dilepton mass consistent 
with the J/ip or ip{2S) mass, the candidate is vetoed. 
There is also significant feed-up from B — > J/ipK and 
B ■)jj{2S)K into B K*l+£-, since energy lost due 
to bremsstrahlung in i? — > J/^K can be compensated 
for by including a random pion. If the Ki~^£~ system 
in a i? ^ K*£^£^ candidate is kinematically consistent 
with B J/^fj{^ £'^£~j)K, assuming that the photon 
(which is not directly observed) was radiated along the 
direction of either lepton, then the candidate is vetoed. 
Apart from the charmonium vetoes, we analyze the full 
mi+i- range. 

Continuum background from non-resonant e+e^ qq 
production is suppressed using a Fisher discriminant [Q, 
a linear combination of the input variables with opti- 
mized coefficients. The variables are R2', cos 0b, the 
cosine of the angle between the B candidate and the 
beam axis in the cm. frame; cos^t, the cosine of the 
angle between the thrust axis of the candidate B meson 
daughter particles and that of the rest of the particles 
in the cm. frame; and rriKi, the invariant mass of the 
if-lepton system, where the lepton is selected according 
to its charge relative to the strangeness of the K^*\ The 
variable mKe helps discriminate against background from 
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FIG. 1: Charmonium veto in the A_B vs. rrii+i- plane for (a) 
B K^'^e^e' and (b) B — > K^*^ fi'^ ii~ . Hatched regions are 
vetoed. The dots correspond to a Monte Carlo simulation of 
B J/ip{^ e+i-)K and B i/'(25')(^ l+r)K. Most sig- 
nal events would lie in the Ai? region between the horizontal 
lines. 



semileptonic D decays, for which niKe. < m^i. 

Combinatorial background from BB events is sup- 
pressed using a signal-to- _B_B hkelihood ratio that com- 
bines candidate B and dilcpton vertex probabilities; the 
significance of the dilepton separation along the beam 
direction; cos 9b] and the missing energy, Smiss, of the 
event in the cm. frame. The variable ii^miss provides the 
strongest discrimination against BB background, since 
events with semileptonic decays usually have significant 
unobserved energy due to neutrinos. For each final state, 
we select at most one combination of particles per event 
as a i? signal candidate. If multiple candidates occur, 
we select the candidate with the greatest number of drift 
chamber and SVT hits on the charged tracks. 

We use the known charmonium decays B — > J/^K^*'^ 
and B — > ^p{2S)K^*'> to check the efficiency of our analy- 
sis cuts. Figure || compares the AE distributions (ab- 
solutely normalized) of these charmonium samples in 
Monte Carlo with data. We find good agreement in both 
the normalization and the shape. 

We extract the signal and background yields in 
each channel using a two-dimensional extended un- 
binned maximum likelihood fit in the region defined by 
TOES > 5.2 GeV/c2 and \AE\ < 0.25 GeV. The signal 
shapes, including the effects of radiation on the AE dis- 
tribution and the correlation between toes and AE, are 
obtained by parametrizing the GEANT3 Monte Carlo ||] 
simulation of the signal. The background is described by 
a function |^ with two parameters that are determined 
in our fits to the data. Backgrounds from BB that peak 
in the signal region are suppressed to less than 0.2 events 
in each mode. Although we allow the signal yield to be 
negative, we have imposed a lower cut-off such that the 
total fit function is positive. The fit results are shown 
in Fig. ^ and summarized in Table || We observe no 
significant signals. 

To determine 90% C.L. upper limits on the signal 



yields, we generate and fit a series of toy Monte Carlo 
samples in which the background probability density 
function is taken from our fit to the data, but the mean 
number of signal events is varied. We generate ten thou- 
sand samples for each mean value, increasing the mean 
until 90% of the fits to a set of samples give a signal yield 
greater than that obtained by fitting the data. To give a 
measure of the sensitivity of the analysis we list in Table || 
an effective background yield. This quantity is defined as 
the square of the error on the signal yield from a fit to 
a toy Monte Carlo sample drawn from the background 
probability function, with no signal contribution. 

Table | lists the systematic uncertainties from the fit, 
{AB/B)fit, expressed according to their effect on the lim- 
its. The sensitivity of the limits to the values used for 
signal-shape parameters is determined by performing al- 
ternative fits using parameters from the B J/^pK^*") 
control samples. For modes with electrons, we also varied 
the fraction of signal events in the tail of the AE distribu- 
tion. To determine whether a more general background 
shape would lead to different results, we introduced addi- 
tional parameters and allowed for a correlation between 
Toes and AE. This procedure shifted the upper lim- 
its by 2% to 5%, depending on the mode. Most of the 
uncertainty associated with the background shape is in- 
corporated in the statistical error on the yield because 
the background shape is determined from the fit. 
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FIG. 2: Comparison of event yields and AE shapes between 
data and Monte Carlo for the charmonium control samples. 
The points with error bars show the data, and the solid his- 
tograms show the prediction of the charmonium Monte Carlo. 
All of the analysis selection criteria have been applied except 
for the charmonium veto, which is reversed. The large tails 
in the e~ modes are due to photon radiation. Small shifts 
between data and Monte Carlo are taken into account as sys- 
tematic uncertainties on the signal yields. 
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TABLE I; Results from the fits to i? — + K^*"^ t^l~ and B A''*'e*/i^ modes. Tfie columns from left to right are fitted signal 
yield upper limit on the signal yield; the contribution of the background to the error on the signal yield, expressed as an 
effective background yield (see text); the signal efficiency, e (not including the branching fractions for K* , K'^ , and Kg decays); 
the systematic error on the selection efficiency, {AB/B)^; the systematic error from the fit, (AB/B)fit; the branching fraction 
central value (B); and the upper limit on the branching fraction, including systematic errors. 





Mode 


Signal 


90% C.L. 


Effective 


e 


(AS/S). iAB/B)m 


B/10-'' 


B/IO-" 






yield 


yield 


background (%) 


(%) 


(%) 




90% C.L. 
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n n+i.b 

'J-^-o.o 


3.1 


0.7 


17.5 


±7.6 


±4.0 
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0.6 


10.5 
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1.2 
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8.8 


1.4 
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U-'J-O.O 


3.5 
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"•^-0.0 


3.2 


B" 
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0.2 


15.7 


±8.8 


±9.5 


n q+2-2 


3.9 


B" 




'-'•"-0.0 
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0.1 


9.6 
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B+ 


^ K*+e+e- 


'-'•^-0.0 
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±5.0 
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3.9tii 
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B+ 
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1.3 


16.8 


±5.7 


±4.0 


-0 1+"-'' 
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±3.0 
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FIG. 3: Projections from individual maximum likelihood ffts 
onto mES for the AE signal regions: —0.11 < AE < 0.05 GeV 
(electrons) and -0.07 < AE < 0.05 GeV (muons). The dot- 
ted lines show the background component, and the solid lines 
show the sum of background and signal components. 



The systematic uncertainties on the efficiency, 
{AB/B)i, are listed in Table || and arise from charged- 
particle tracking (±1.2%/lepton, ±2.0% for the pion from 
K* Ktt, and ±1.3%/track for other charged hadrons), 
particle identification (±1.4%/electron, ±1.0%/muon, 
±2.0%/track for kaons and pions), the continuum sup- 



pression cut (±2.0%), the BB suppression cut (±3.0%), 
Kg selection (±4.0%), Monte Carlo signal statistics 
(±3.0% to ±5.0%), the theoretical model dependence of 
the efficiency (±4.0% to ±7.0%, depending on the mode), 
and the number of BB events (±1.6%). The uncertain- 
ties on the efficiencies due to mo del- dependence of form 
factors are taken to be the full range of variation ob- 
tained from different theoretical models In setting 
an upper limit, the systematic uncertainties from the ef- 
ficiency, {AB/B)e, and from the fit, {AB/B)fit, are added 
in quadrature, and the limit is increased by this factor. 

Table || also includes the results for the lepton-family- 
violating decays B — > K^*'>efj,, where the signal efficien- 
cies were determined from phase-space Monte Carlo sim- 
ulations. We observe no evidence for these decays. 

We determine the branching fractions B{B ^ Ki'^l^) 
and B{B K*£^£~) averaged over both B meson 
charge and lepton type (e^e~ and fX^ fjT) by perform- 
ing a simultaneous maximum likelihood fit to the four 
contributing channels in each case. In combining the 
B K*£^£^ modes, the ratio of branching fractions 
B{B ^ K*e+e-)/B{B K*^i+fi-) = 1.2 from the 
model of Ali et al. [Q is used to weight the yield in the 
muon channel relative to that in the electron channel. 
The extracted yield corresponds to the electron mode. 
The combined fits give 



B{B 
B{B- 



-> K£+t 
K*£+£- 



= (-0.06l|]:oo ± 0.03) X 10"^ 



(0.9: 



-1.3 
-0.9 



-0.00 

±0.1) X 10" 



where the first error is statistical and the second is sys- 
tematic. We evaluate the upper limits on these combined 
modes and obtain 



B(B -> K£+i- 
B{B ^ K*£+t^ 



< 0.50 X 10-^ at 90% C.L. 

< 2.9 X 10"^ at 90% C.L. 
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These limits represent an improvement over previously 
published results from CDF and CLEO The 
Belle experiment has also recently obtained results 
on these modes. We see no evidence for a signal, and 
our limits are close to many of the predictions based on 
the Standard Model. With the rapidly increasing size of 
our data sample, we expect to have significantly better 
sensitivity to these modes in the future. 
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